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Abstract. We discuss vector boson production at hadron colliders and the treatment 
J> \ of the initial-state radiation according to Monte Carlo parton shower simulations and 

resummed calculations. In particular, we investigate the effect of matrix-element cor- 
| rections to the HERWIG event generator on W/Z transverse momentum distributions. 
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The production of vector bosons W, Z and 7 [1] is one of the most interesting 
processes in the phenomenology of hadron collisions and provides an environment 
to test both Quantum Chromodynamics and the Standard Model of electroweak 
interactions (see [2,3] for a review). The lowest-order processes qq' — > V are not 
sufficient to make reliable predictions, but the initial-state radiation should be taken 
into account. Monte Carlo event generators and resummed analytical calculations 
are available tools to describe the multiple radiation accompanying the incoming 
hadrons. 

Standard Monte Carlo algorithms [4,5] describe the initial-state parton showers 
in the soft /collinear approximation, but can have 'dead zones', where no radiation is 
permitted. The radiation in these regions can be generated by the use of the exact 
first-order matrix element. Referring to the HERWIG event generator, matrix- 
element corrections to Drell-Yan processes have been implemented in [6] , following 
the general method of [7], and included in the new version HERWIG 6.1 [8]. 

Another possible approach consists of performing an analytical resummation of 
the large logarithmic coefficients which multiply the strong coupling constant. Con- 
sidering the transverse momentum q? distribution, logarithms of the ratio rriv/qT, 
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m v being the vector boson mass, arise in calculating higher-order corrections to 
the Born process. The resummation of these logarithms, which are large at small 
qT, was initially proposed by Dokshitzer, Dyakonov and Troyan (DDT) [9], then 
accomplished by Collins, Sterman and Soper (CSS) [10]. CSS performed the resum- 
mation in the space of the impact parameter b, which is the Fourier conjugate of 
qx- Their results have been implemented numerically in [11,12], while more recent 
analyses can be found in [13-15], where the resummation is performed in both q T - 
and 6-space. 

In this paper, we review some results for the W/Z transverse momentum distri- 
bution according to the HERWIG event generator and resummed calculations. 2 



HERWIG simulates the initial-state radiation in hadron collisions according to 
a 'backward evolution' [17], in which the scale is reduced away from the hard ver- 
tex and traces the hard-scattering partons back into the incoming hadrons. The 
branching algorithm relies on the universal structure of the elementary probabil- 
ity in the leading infrared approximation. The probability of the emission of an 
additional soft /collinear parton from a parton i is given by the general result: 



and £j = with p and pi being the four- momenta of the splitting and of the 
emitted parton respectively; Z{ is the energy fraction of the outgoing space-like 
parton with respect to the incoming one; P a b{z) is the Altarelli-Parisi splitting 
function for a parton a evolving in b. In the approximation of massless partons, 
we have = 1 — cos#, where is the emission angle to the incoming hadron 
direction. For soft emission (Ei <C E), ordering according to qf corresponds to 
angular ordering. When the emission is hard, the energy of the radiated parton 
is similar to that of the splitting parton, so qf-ordering corresponds to transverse 
momentum ordering. In (1) f a (xi,qf) is the parton distribution function for the 
partons of type a in the initial-state hadron, being the parton energy fraction. 
The function 



is the Sudakov form factor, expressing the probability of no-resolvable branching in 
the range q\ < q 2 < q\. The ratio of form factors in (1) is therefore the probability 
of no branching at higher values of qf. Unitarity dictates that the Sudakov form 
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factor sums up all-order virtual and unresolved contributions. In (1), gj max is the 
maximum value of q, fixed by the hard process, and q c is the value at which the 
backward evolution is terminated, corresponding, in the case of HERWIG, to a 
cutoff on the transverse momentum of the showering partons. However, since q c 
is smaller than the minimum scale at which the parton distribution functions are 
evaluated, an additional cutoff on the evolution variable q 2 has to be set. 

If the backward evolution has not resulted in a valence quark, an additional non- 
perturbative parton emission is generated to evolve back to a valence quark. Such a 
valence quark has a Gaussian distribution with respect to the non-perturbative in- 
trinsic transverse momentum in the hadron, with a width q Tint that is an adjustable 
parameter and whose default value is zero. 

We need finally to specify the showering frame, the variables q 2 and Zi being 
frame-dependent. One can show that, as a result of the g 2 -ordering, the maximum 
q- values of two colour connected partons % and j are related Via Q'imaxQ'imax Pi 'Pji 
which is Lorentz- invariant. For vector boson production, symmetric limits are set 
in HERWIG: qi max = qj max = y/Pi • Vj ■ Furthermore, the energy of the parton 
which initiates the cascade is given by E = q max = v /pi _r pj. It follows that ordering 
according to q 2 implies £ < z 2 . 

The region £ > z 2 is therefore a 'dead zone' for the shower evolution. In such 
a zone the physical radiation is not logarithmically enhanced, but not completely 
absent as happens in the standard algorithm. We therefore need to improve the 
HERWIG parton showers by the use of matrix-element corrections. 

MATRIX-ELEMENT CORRECTIONS 

According to [7], we populate the 'dead zone' of the phase space using the ex- 
act 0(ots) matrix element (hard correction). We also correct the emission in the 
already-populated region using the first-order result any time an emission is capa- 
ble of being the 'hardest so far' (soft correction), where the hardness of an emission 
is measured in terms of the transverse momentum of the emitted parton relative 
to the splitting one. 

We consider the process q{pi)q'{p2) V(q)9(P3)j define the Mandelstam vari- 
ables s — (pi + P2) 2 , t — (pi — ps) 2 and u = (p 2 — p%) 2 and obtain the total 
phase-space limits 

m v < s< s, (3) 
my - s < i< 0, (4) 

s being the total centre-of-mass energy. We observe that the soft singularity corre- 
sponds to s = my and the lines t — and t = m v — s to collinear gluon emission. 

After relating the parton shower variables z and £ to s and t, as done in [6], and 
setting £ < z 2 , one can get the HERWIG phase-space limits in terms of s and i. 
In Fig. 1 we plot the total and the HERWIG phase space for y/s = 200 GeV and 
my = 80 GeV; the soft and collinear singularity are inside the HERWIG region. 




FIGURE 1. Total (solid line) and HERWIG (dashed) phase space limits for = 200 GeV and 
m v = 80 GeV. 



We calculate the differential cross section with respect to s and t following the 
prescription of [18], where it is shown that, assuming that the rapidity and the vir- 
tuality of the vector boson are fixed by the Born process, the following factorization 
formula holds: 

,2 fq/l(Xl) fq' /2(X2) C F a S d,S dt r. 2 2 2 ; 2 1 /r\ 

d a = a ° t 7Z, VP („ \ ^ ( m v- u ) +{m v -t) , (5) 

J q /l{Vl> Jq' /2{V2> 27T SHU L J 

where f q /i(xi) an d fq> /2(X2) are the parton distribution functions of the scattering 
partons inside the incoming hadrons 1 and 2 for energy fractions Xi an d X2 in the 
process qq' — > while f q /i(r)i) and fqi^fa) refer to the Born process. 

A similar treatment holds for the Compton process qij>i)g(j>z) Q f (P2)V (q) , as 
discussed in [6]. 

The distribution (5) or the equivalent one for the Compton process is imple- 
mented to generate events in the missing phase space and in the populated region 
every time an emission is the hardest so far. 



TRANSVERSE MOMENTUM DISTRIBUTIONS: 
HERWIG RESULTS 

An interesting observable to study is the vector boson transverse momentum, 
which is constrained to be q T < rn v in the soft /collinear approximation. After 
matrix-element corrections, a fraction of events at higher q^ is to be expected. In 
Figs. 2 and 3 we plot the W qx distribution at the Tevatron and at the LHC, 
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FIGURE 2. W transverse momentum distributions at the Tevatron, according to HERWIG 6.1 
(solid line) and HERWIG 5.9 (dotted). 

according to HERWIG 5.9 and HERWIG 6.1, the new version including matrix- 
element corrections, for q Tint = 0. We see a big effect at large qx- after some qx the 
5.9 version does not generate events anymore, while we still have a non-zero cross 
section after matrix-element corrections. 

Moreover, a slight suppression can be seen at small qx- It is related to the fact 
that, although we are providing the shower with the tree-level 0(as) corrections, 
virtual contributions are missing and we still get the leading-order cross section. 
The enhancement at large q T is therefore compensated by a suppression in the 
low-g-r range. 

It is now interesting to compare the HERWIG results with some Tevatron data. 
In Fig. 4 we compare the HERWIG 6.1 distribution with some D0 data [19] and 
find reasonable agreement over the whole qx range. As shown in [6], smearing the 
HERWIG curve to account for detector effects must be included to achieve this 
agreement. Also, we do not see any relevant impact of setting qxmt — 1 GeV after 
detector corrections. In Fig. 5 we compare the HERWIG 5.9 and 6.1 results, for 
different values of (/Tint, with some CDF data [20], already corrected for detector 
effects. We find good agreement after matrix-element corrections, while the 5.9 
version is not able to fit in with the data for qx > 50 GeV. At low qx, the best 
agreement to the data is obtained for qx\nt = 2 GeV, as shown in Fig. 6. While 
the Z distribution strongly depends on g Tint at small qx, in [21] and Fig. 7 it is 
shown that the ratio R of the W and Z spectra is approximately independent of 
it. 3 . Such a ratio is one of the main inputs for the W mass measurement in hadron 
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FIGURE 3. As in Fig. 2, but at the LHC. 




FIGURE 4. Comparison of the D0 data to the HERWIG 6.1 result, for an intrinsic q T of 
(solid line) and 1 GeV (dashed). 
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FIGURE 5. Comparison of the CDF data on the Z transverse momentum to HERWIG 5.9 
(dotted line) and 6.1 for an intrinsic qr of (solid line), 1 GeV (dashed) and 2 GeV (dot-dashed). 

collisions and it is good news that it does not depend on unknown non-perturbative 
parameters. 



Another possible approach to study the vector boson transverse momentum dis- 
tribution consists of resumming the logarithmic terms / = \og(m v /q T ) in the low- 
qT range. It is interesting to compare the HERWIG phenomenological results with 
those of some resummed calculations, in particular [13] and [14]. According to [9], 
the resummed differential cross section for W production can be written as: 



where V qq -> is the relevant Cabibbo-Kobayashi-Maskawa matrix element and r = 
rriy/s. In (6), exp[5 , (my, qr)} is a Sudakov-like form factor which resums the large 
logarithms associated to the initial-state radiation. It reads: 



RESUMMED CALCULATIONS 



drriydq; 




(6) 



S(m v , qr) = - U s (/i 2 )) log + B(a s ^ 2 )) 



(7) 



where A{ag) and B(pts) can be expanded as: 





FIGURE 7. The ratio of the W to the Z spectrum at small qx for an intrinsic transverse 
momentum of zero (solid), 1 (dashed) and 2 GeV (dotted). 



A(as) = Axas + A 2 a 2 s + . . . ; B(a s ) = B ± as + B 2 a 2 s + 



(8) 



As far as the logarithms which contribute to the resummation are concerned, two 
conflicting nomenclatures exist. One consists of looking at Sudakov exponent, 
where the leading logarithms (LL) are ~ agl n+1 and the next-to-leading ones (NLL) 
~ oJgl 71 . It is straightforward to show that the LL contributions are obtained by 
keeping only the A\ term in the expansions (8) while NLL accuracy is achieved by 
considering A 2 and B\ as well. In this sense, the approach [14] is NLL. 
Another classification relies on the expansion of the exponent 

S(m v , q T ) = J2 c n ,n + ia n s l n+1 + £ c n , n a n s l n , (9) 

n n 

where the leading term is ~ cts^ 2 , so that the leading contributions to 
exp[S (my , qr)} are ~ agl 2n , terms ~ a^l 2 " 1 ^ 1 being next-to-leading. This is equiv- 
alent to saying that in the differential cross section the LL and NLL contributions 
are ~ (1/ 'q T )agl 2n ~ 1 and ~ (l/q' T )agl 2n ~ 2 respectively. According to this nomen- 
clature, the calculations [13] and [15] are NNLL and NNNLL respectively. 
In the 6-space formalism, following [14], the differential cross section reads: 



d 2 cr a 



2 iKd 2 f dx 1 dx 2 S( Xl x 2 -t) J d 2 b e i&X 



dm\dq\ 4vr ^ Jo 

x [f g /i(xi, c 1 /b)fq //2 (x 2 , Ci/6) exp[S(m v , b)} + (q <-> q )] , (10) 

where c\ and c 2 are integration constants of order 1 and S(m v ,b) is the Sudakov 
exponent in 6-space. 

For high b values, i.e. small qT, non-perturbative effects are taken into account 
via a Gaussian function F^p = exp(—gb 2 ), as suggested in [11]. In both [13] and 
[14] the value g = 3 GeV 2 is chosen. 

Also, in order to allow resummed calculations to be reliable even at large qr, we 
wish to match the calculations of [13] and [14] to the exact 0(ocs) result. We add 
the first-order cross section to the resummed result and, in order to avoid double 
counting, we subtract off the term which they have in common, which is the qx — > 
limit of the exact 0{ag) result. According to our prescription, the matching works 
fine if at the point qr = my we have a continuous distribution. 



COMPARISON OF HERWIG AND RESUMMED 

CALCULATIONS 

A detailed and general discussion on the comparison of angular-ordered parton 
shower algorithms with resummed calculations for Drell-Yan processes was already 
performed in [22], where the authors showed that, for r — > 1, HERWIG always 
accounts for the term Ai, corresponding to the leading logarithms in the exponent, 
and B\ as well. Furthermore, one is able to account for the NLL term A 2 by 



simply modifying the Altarelli-Parisi splitting function introducing a second-order 
contribution 



P' qq {a s ,z) = ^C F - 
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where the K factor is given by: 
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Nf being the number of flavours, C F = 4/3 and Ca = 1/2. This is equivalent to 
redefining the QCD parameter A to the 'Monte Carlo' A M c- 



with p = (HC A - 2JV / )/(127r). Even after these replacements, the HERWIG 
algorithm cannot be considered completely accurate at the next-to-leading level, 
since it is still missing higher-order contributions in the strong coupling constant 
or the parton distribution functions (see, for instance, the discussion in [3]). 

In Fig. 8 we show the W transverse momentum distribution at the Tevatron 
in the low-qx range according to HERWIG 6.1 and the calculations of [13] in qx- 
space and of [14] in q T - and 6-space. The HERWIG curve lies within the range of 
the resummed calculations, which is a reasonable result, considering that we are 
actually comparing different approaches. In Fig. 9 we consider the whole qx range, 
with the resummations matched to the exact first-order amplitude. We find that 
the matching works fine only for the approach [14] in g T -space, the others showing 
a step at qx = m\v, due to uncompensated contributions of order a 2 s or higher. 
The well-matched distribution agrees with the HERWIG 6.1 prediction at large qx- 



We studied the initial-state radiation in vector boson production according to the 
HERWIG event generator and some resummed calculations. In particular, we in- 
vestigated the effect of the recently-implemented matrix-element corrections to the 
HERWIG algorithm. We found a big effect of such corrections on W/Z transverse 
momentum distributions at the Tevatron and at the LHC, and good agreement 
with the D0 and CDF data, with a crucial role played by such corrections in order 
to be able to fit in with the data at large qx- We also found that, even though 
the spectra at small qx do depend on the intrinsic non-perturbative transverse 
momentum, the ratio of the W to the Z spectrum is roughly independent of it. 
We then considered some resummed calculations, which we matched to the exact 
0(as) matrix element, which makes them reliable at large qx as well. We found 
reasonable agreement of such approaches with HERWIG and fine matching only 



Kmc = Aexp(A747rA)), 



(13) 



CONCLUSIONS 
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FIGURE 8. The W qr distribution at the Tevatron, according to HERWIG with matrix- element 
corrections, with zero intrinsic qx (solid histogram) and an qnnt of 1 GeV (dashed histogram), 
compared with the resummed results of [14] in g^-space (solid line) and in &-space (dotted line) 
and of [13] (dashed line) in qx-sp&ce. 
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FIGURE 9. As in Fig. 8, but with the resummed results matched with the exact 0(as) matrix 
element. 



for the calculation which keeps all the next-to-leading logarithms in the Sudakov 
exponent in g^-space. 

Finally, we have to say that the discussed method of improving the initial-state 
radiation in parton-shower Monte Carlo simulations can be extended to a wide 
range of interesting processes for the phenomenology of hadron colliders. The 
implementation of hard and soft corrections to top and Higgs production is in 
progress. 
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